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Abstract 

The third Newton law for mean velocity fields is utilised to gener- 
ate anomalous (enhanced) or non-dispersive diffusion-type processes 
which, in particular, can be interpreted as a probabilistic counterpart 
of the Schrodinger picture quantum dynamics. 



If we consider a fluid in thermal equilibrium as the noise carrier, a kinetic 
theory viewpoint amounts to visualizing the constituent molecules that col- 
lide not only with each other but also with the tagged (colloidal) particle, so 
enforcing and maintaining its observed erratic motion. The Smoluchowski 
approximation takes us away from those kinetic theory intuitions by project- 
ing the phase-space theory of random motions into its configuration space 
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image which is a spatial Markovian diffusion process, whose formal infinites- 
imal encoding reads: 

dX(t) = -^=dt + V2DdW(t) . (1) 

mp 

In the above m stands for the mass of a diffusing particle, (3 is a friction 
parameter, D is a diffusion constant and W(t) is a normalised Wiener process. 
The Smoluchowski forward drift can be traced back to a presumed selective 
action of the external force F = — VV on the Brownian particle that has 
a negligible effect on the thermal bath but in view of frictional resistance 
imparts to a particle the mean velocity F/m/3 on the (3~ l time scale, [p]]. 
The noise carrier (fluid in the present considerations) statistically remains in 
the state of rest, with no intrinsic mean flows. 

An implicit phase-space scenario of the Brownian motion refers to minute 
acceleration/deceleration events which modify (infinitely gently, [|lj, at a 
generic rate of 10 21 times per second) velocities of realistic particles. Clearly, 
the microscopic energy-momentum conservation laws need to be respected 
in each separate collision event. In contrast to derivations based on the 
Boltzmann equation, this feature is completely alien to the Brownian mo- 
tion theory. That energy-momentum deficit is one of "forgotten" or tacitly 
"overlooked" problems in the standard theory of the Brownian motion, see 
e.g. §!■ 

An issue of origins of general isothermal flows in the noisy environment is 
far from being settled, || . Namely, the standard theory does not account (or 
refers to phenomena for which that is not necessary) for generic perturbations 
of the random medium as a reaction to the enforced (exclusively due to the 
impacts coming from the medium) particle motion. An ability of the medium 
to perform work i. e. to give a kinetic energy to Brownian particles (strictly 
speaking, in the local mean, hence on the ensemble average) appears to be 
a universal feature of the thermal bath even in the absence of any external 
force. 
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Then, there arise problems with the thermal equilibrium, ||. The tagged 
particle always propagates "at the expense" of the bath, which nonetheless 
remains " close" to its thermal equilibrium. In view of the involved local cool- 
ing and heating phanomena, which are necessary to keep in conformity with 
the second law of thermodynamics, the bath should develop local flows and 
thus in turn actively react back to what is being happenning to the particle 
in the course of its random propagation. Such environmental reaction, while 
interpreted as a feedback effect, surely can be neglected in an individual 
particle propagation sample. However, an obvious Brownian particle energy- 
momentum "deficit" on the ensemble average (kinetic energy and momentum 
is carried locally by diffusion currents), indicates that the medium reaction 
may give observable contributions on the ensemble average, as a statistically 
accumulated feature. For a statistical ensemble of weakly-out-of-equilibrium 
systems, an effective isothermal scenario can be used again, but then mean 
flows in the noise carrier are "born" as a consequence of the averaging, ||. 

It is well known that a spatial diffusion (Smoluchowski) approximation 
of the phase-space process, allows to reduce the number of independent lo- 
cal conservation laws (cf. to two only. Therefore the Fokker-Planck 
equation can always be supplemented by another (independent) partial dif- 
ferential equation to form a closed system. 

If we assign a probability density po(x) with which the initial data xq = 
X(0) for Eq. (1) are distributed, then the emergent Fick law would reveal a 
statistical tendency of particles to flow away from higher probability residence 
areas. This feature is encoded in the corresponding Fokker-Planck equation: 



flow that might be experimentally observed for a cloud of suspended particles 



d tP = -V • (vp) = -V • [(^- - D^)p] (2) 
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in a liquid) j = vp gives rise to a non- negligible matter transport on the 
ensemble average, [|]. ||]. 

It is interesting to notice that the local velocity field v(x, t) obeys the 
natural (local) momentum conservation law which directly originates from 
the rules of the ltd calculus for Markovian diffusion processes, and from 
the first moment equation in the diffusion approximation (!) of the Kramers 
theory, [|: 

d t v + (v- V)if = V(fi-Q) . (3) 

An effective potential function Q(x) can be expressed in terms of the 
forward drift b(x) = as follows: f2 = „ F l n2 + -^V • F. 

\ / mp 2m^p z mp 

Let us emphasize that it is the diffusion (Smoluchowski) approximation 
which makes the right- hand- side of Eq. (3) substantially different from the 
usual moment equations appropriate for the Brownian motion. In particular, 
the force F presumed to act upon an individual particle, does not give rise 
in Eq. (3) to the expression — — W which might be expected on the basis of 
kinetic theory intuitions and moment identities directly derivable from the 
Karmers equation, but to the term +Vfi. 

Moreover, instead of the standard pressure term, there appears a contri- 
bution from a probability density p-dependent potential Q(x,t). It is given 
in terms of the so-called osmotic velocity field u(x,t) = D'V lnp(x,t), (cf. 
[0]) : Q{%> t) = \u 2 + .DV • u and is generic to a local momentum conservation 
law respected by isothermal Markovian diffusion processes, cf. [fjj-|§. 

The Smoluchowski drift does not refer to any flows in the noise carrier. To 
analyze perturbations of the medium and then the resulting intrinsic (mean) 
flows, a more general function b(X(t),t), must replace the Smoluchowski drift 
in Eqs. (1), (2). Forward drifts modify additively the pure noise (Wiener 
process entry) term in the Ito equations. Under suitable restrictions, we can 
relate probability measures corresponding to different (in terms of forward 
drifts !) Fokker-Planck equations and processes by means of the Cameron- 
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Martin-Girsanov theory of measure transformations. The Radon-Nikodym 
derivative of measures is here involved and for suitable forward drifts which 
are are gradient fields that yields, M, the most general form of an auxiliary 
potential Q(x,t) in Eq. (3): 

Q(x,t)=2D[d t <f ) +^ + V-b)} . (4) 

We denote b(x,t) = 2DW(j)(x,t). 

Eq. (4) is a trivial identity, if we take for granted that all drifts are known 
from the beginning, like in case of typical Smoluchowski diffusions where the 
external force F is a priori postulated. We can proceed otherwise and, on 
the contrary, one can depart from a suitably chosen space-time dependent 
function Q(x,t). From this point of view, while developing the formalism, 
one should decide what is a quantity of a primary physical interest: the field 

— * 

of drifts b(x,t) or the potential Q(x,t). 

Mathematical features of the formalism appear to depend crucially on 
the properties (like continuity, local and global boundedness, Rellich class) 
of the potential Q, see e.g. (|, ||. Let us consider a bounded from below 
(local boundedness from above is useful as well), continuous function fl(x, t), 
cf. [|J). Then, by means of the gradient field ansatz for the diffusion current 
velocity (v = VS — > d t p = — V ■ [(VS')p]) we can transform the momen- 
tum conservation law (3) of a Markovian diffusion process to the universal 
Hamilton- Jacobi form: 

n = d t S + ^\VS\ 2 + Q (5) 

where Q(x,t) was defined before. By applying the gradient operation to 
Eq. (5) we recover (3). In the above, the contribution due to Q is a direct 
consequence of an initial probability measure choice for the diffusion process, 
while Q via Eq. (4) does account for an appropriate forward drift of the 
process. 
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Thus, in the context of Markovian diffusion processes, we can consider a 
closed system of partial differential equations which comprises the continuity 
equation d t p = — V(-up) and the Hamilton- Jacobi equation (5), plus suitable 
initial (and/or boundary) data. The underlying diffusion process is specified 
uniquely, cf. [|, ||. 

Since the pertinent nonlinearly coupled system equations looks discour- 
aging, it is useful to mention that a linearisation of this problem is provided 
by a time-adjoint pair of generalised diffusion equations in the framework of 
the so-called Schrodinger boundary data problem, ||, |^| . The standard heat 
equation appears as a very special case in this formalism. 

The local conservation law (3) acquires a direct physical meaning (the 
rate of change of momentum carried by a locally co-moving with the flow 
volume, only if averaged with respect to p(x,t) over a simply connected 
spatial area. If V stands for a volume enclosed by a two-dimensional out- 
ward oriented surface dV, we define a co-moving volume on small time scales, 
by deforming the boundary surface in accordance with the local current ve- 
locity field values. Then , let us consider at time t the displacement of 
the boundary surface dV(t) defined as follows: x G dV — > x + v(x,t)At 
for all x G dV. Up to the first order in At this guarantees the con- 
servation of mass (probability measure) contained in V at time t i. e. 
Iv(t+At) P& t + &t)d 3 x - J v(t) p(x, t)d 3 x ~ 0. 

The corresponding (to the leading order in At) quantitative momentum 
rate-of-change measure reads, cf. ||J], J v pV(fi — Q)a?x. 

For a particular case of the free Brownian expansion of an initially given 
Po(x) = (TO 2)3/2 egp(-^), where a 2 = ADt , we would have - J v pVQd 3 x = 
— J dv Pda, where Q(x,t) = gT^pp — 2(t+t ) ' wm ^ e ^ e "osmotic pressure" 
contribution reads P(x, t) = — 2 (t+t ) P$i ^ or a ^ % ^ R 3 an d t > 0. 

The current velocity v(x, t) = VS(x, t) = 2 (t+t ) * s nm< ed to the Hamilton- 
Jacobi equation d t S + llVS 1 ) 2 + Q = whose solution reads: S(x,t) = 
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3 ^ + |DZ„[47rI>(t + to)]. 

Let us observe that the initial data vq = —DVlnpo = —uq for the cur- 
rent velocity field indicate that we have totally ignored a crucial preliminary 
stage of the dynamics on the time scale, when the Brownian expansion 
of an initially static ensemble has been ignited and so particles have been 
ultimately set in motion. 

Notice also that our "osmotic expansion pressure" P(x,t) is not posi- 
tive definite, in contrast to the familiar kinetic theory (equation of state) 
expression for the pressure P(x) = q /(i), a > appropriate for gases. 
The admissibility of the negative sign of the "pressure" function encodes 
the fact that the Brownian evolving concentration of particles generically 
decompresses (blows up), instead of being compressed by the surrounding 
medium. 

The loss (in view of the "osmotic" migration) of momentum stored in 
a control volume at a given time, may be here interpreted in terms of an 
acceleration — J v pVQd 3 x induced by a fictituous " attractive force" . 

By invoking an explicit Hamilton- Jacobi connection (5), we may attribute 
to a diffusing Brownian ensemble the mean kinetic energy per unit of mass 
J v p\v 2 d 3 x. In view of < x 2 >= 6D(t + to)) we have also f R 3 p\v 1 d 3 x = 
4 u+ to y Notice that the mean energy f v p(^v 2 + Q)d 3 x needs not to be posi- 
tive. Indeed, this expression identically vanishes after extending integrations 
from V to R 3 . On the other hand the kinetic contribution, initially equal 
J R 3 \pv 2 d 3 x = 3D /a 2 and evidently coming from nowhere, continually di- 
minishes and is bound to disappear in the asymptotic t — > oo limit, when 
Brownian particles become uniformly distributed in space. 

Normally, diffusion processes yielding a nontrivial matter transport (dif- 
fusion currents) are observed for a non-uniform concentration of colloidal par- 
ticles which are regarded as independent (non- interacting). We can however 
devise a thought (numerical) experiment that gives rise to a corresponding 
transport in terms of an ensemble of sample (and thus independent) Brown- 
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ian motion realisations on a fixed finite time interval, instead of considering a 
multitude of them (migrating swarm of Brownian particles) simultaneously. 

Let us assume that "an effort" (hence, an energy loss/gain which implies 
local deviations from thermal equilibrium conditions) of the random medium, 
on the f3~ l scale, to produce a local Brownian diffusion current (pv)(x,t ) 
out of the initially static ensemble and thus to decompress (lower the blow- 
up tendency) an initial non-uniform probability distribution, results in the 
effective osmotic reaction of the random medium. This is the Brownian recoil 
effect of Ref. [|. 

In that case, the particle swarm propagation scenario becomes entirely 
different from the standard Brownian one. First of all, the nonvanishing 
forward drift b = u is generated as a dynamical (effective, statistical here !) 
response of the bath to the enforced by the bath particle transport with the 
local velocity v = —u. Second, we need to account for a parellel inversion 
of the pressure effects (compression +VQ should replace the decompression 

— VQ) in the respective local momentum conservation law. 

Those features can be secured through an explicit realization of the action- 
reaction principle which we promote to the status of the third Newton law in 
the mean. 

On the level of Eq. (3), once averaged over a finite volume, we interpret 
the momentum per unit of mass rate-of-change f v pV(Q—Q)d 3 x which occurs 
exclusively due to the Brownian expansion mechanism, to generate a coun- 
terbalancing rate-of-change tendency in the random medium. To account 
for the emerging forward drift and an obvious modification of the subse- 
quent expansion of an ensemble of particles, we re-define Eq. (3) by setting 

— J v pV(Q — Q)d 3 x in its right-hand-side instead of + J v pV(Q — Q)d 3 x . 
That amounts to an instantaneous implementation of the third Newton law 
in the mean (action-reaction principle) in Eq. (3). 

Hence, the momentum conservation law for the process with a recoil 
(where the reaction term replaces the decompressive "action" term) would 
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read: 



d t v + (v ■ V)v 



V(Q - n) 



(6) 



so that 




(7) 



stands for the corresponding Hamilton- Jacobi equation, cf. [^], [TTJ, instead 
of Eq. (5). A suitable adjustment (re-setting) of the initial data is here 
necessary. 

In the coarse-grained picture of motion we thus deal with a sequence of 
repeatable scenarios realised on the Smoluchowski process time scale At: the 
Brownian swarm expansion build-up is accompanied by the parallel counter- 
flow build-up, which in turn modifies the subsequent stage of the Brownian 
swarm migration (being interpreted to modify the forward drift of the pro- 
cess) and the corresponding built-up anew counterflow. 

The new closed system of partial differential equations refers to Markovian 
diffusion-type processes again, |], £|, [|. The link is particularly obvious if we 
observe that the new Hamilton- Jacobi equation (7) can be formally rewritten 
in the previous form (5) by introducing: 



where f2 r = 2Q—VL and Vt represents the previously defined potential function 
of any Smoluchowski (or more general) diffusion process. 

It is Q r which via Eq. (4) would determine forward drifts of the Marko- 
vian diffusion process with a recoil. They must obey the Cameron-Martin- 
Girsanov identity Q r = 2Q - Q = 2D[d t (f> + |(J^ + V • &)]. 

Our new closed system of equations is badly nonlinear and coupled, but 
its linearisation can be immediately given in terms of an adjoint pair of 
Schrodinger equations with a potential Q, [[[], Indeed, id t ip = —DAip + 
■^ip with a solution ip = p l l 2 exp(iS) and its complex adjoint makes the job, 
if we regard p together with 5* to remain in conformity with the previous 



n r = d t S+-\VS\ 2 + Q 



(8) 
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notations. The choice of ip{x, 0) gives rise to a solvable Cauchy problem. 
(Notice that by setting D = we recover the standard quantum mechanical 
notation.) 

This feature we shall exploit in below. Notice that, for time-indepedent Q, 
the total energy /^(y — Q + Q)pd 3 x of the diffusing ensemble is a conserved 
quantity. 

The general existence criterions for Markovian diffusion processes of that 
kind, were formulated in Ref. 0, see also |5], ff|. 

Let us consider a simple one-dimensional example. In the absence of ex- 
ternal forces, we solve the equations (in space dimension one) d t p = —V(vp) 
and d t v + {vV)v = +VQ, with an initial probability density po(x) cho- 
sen in correspondence with the previous free Brownian motion example. 
We denote a 2 = 4Dt - Then, p(x,t) = [7T{a4+4 a D 2 t2)] i/2 exp[- J^&t* ] and 
b(x,t) = v(x,t) + u(x,t) = 2D }a + ^2i2 X are the pertinent solutions. Notice 
that u(x, 0) = — 2 j^f = b(x, 0) amounts to v(x, 0) = 0, while in the previous 
free Brownian case the initial current velocity was equal to —DVln p . This 
re-adjustment of the initial data can be interpreted in terms of the counterbal- 
ancing (recoil) phenomenon: the would-be initial Brownian ensemble current 
velocity Vq = —Uq is here completely saturated by the emerging forward drift 
b = Uq, see e.g. also @. We deal also with a fictituous "repulsive" force, 
which corresponds to the compression (pressure upon) of the Brownian en- 
semble due to the counter-reaction of the surrounding medium. We can write 
things more explicitly. Namely, now: Q(x,t) = J^l t2 ( ^qifca -1) and the 
corresponding pressure term (VQ = ^VP) reads P(x,t) = — a 4^^t t 2 p(x, t) 
giving a positive contribution +VQ to the local conservation law. The re- 
lated Hamilton- Jacobi equation d t S + ^\VS\ 2 = +Q is solved by S(x,t) = 
Ji+lpiti — D arctan (— With the above form of Q(x, t) one can readily 
check that the Cameron-Martin-Girsanov constraint euqation for the forward 
drift of the Markovian diffusion process with a recoil is automatically valid 
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for = \lnp + S: 2Q = 2D[d t (f> + \{^ + V • 6)]. 

In anology with our free Brownian motion discussion, let us observe that 
presently < x 2 >= ^ + 2D ^f . It is easy to demonstrate that the quadratic 
dependence on time persists for arbitrarily shaped initial choices of the prob- 
ability distribution po(x) > 0. That signalizes an anomalous behaviour (en- 
hanced diffusion) of the pertinent Markovian process when Q = i. e. 

— > 

We can evaluate the kinetic energy contribution § R p^dx = ^T^j^F) 
which in contrast to the Brownian case shows up a continual growth up to 
the terminal (asymptotic) value This value was in turn an initial kinetic 
contribution in the previous free Brownian expansion example. In contrast 
to that case, the total energy integral is now finite (finite energy diffusions 
of Ref. 0) and reads S R (\v 2 — Q)pdx = (it is a conservation law). The 
asymptotic value of the current velocity v ~ f is twice larger than this 
appropriate for the Brownian motion, v ~ 

It is easy to produce an explicit solution to (7), (8) in case of Q(x) = 
}pl 2 x 1 — -D7, with exactly the same inital probability density po(x) as before. 
The forward drift of the corresponding diffusion-type process does not show 
up any obvious contribution from the harmonic Smoluchowski force. It is 
completely eliminated by the Brownian recoil scenario. One may check that 
b(x, 0) = — "^f = u(x,0), while obviously b = ^ = — 7X would hold true 
for all times, in case of the Smoluchowski diffusion process. Because of the 
harmonic attraction and suitable initial probability measure choice, we have 
here wiped out all previously discussed enhanced diffusion features. Now, 
the dispersion is attentuated and actually the non-dispersive diffusion-type 
process is realised: < x 2 > does not spread at all despite of the intrinsically 
stochastic nature of the dynamics (finite-energy diffusions of Ref. 0). 

It is clear that stationary processes are the same both in case of the 
standard Brownian motion and the Brownian motion with a recoil. The 
respective propagation scenarios substantially differ in the non-stationary 



11 



case only. 

One may possibly argue, that a sign inversion of the right-hand-side of 
Eq. (3) which implies (6), can be accomplished by means of an analytic 
continuation in time (cf. the Euclidean quantum mechanics discussion in 
Ref. H). An examination of free dynamics cases given in the above proves 
that the third Newton law ansatz is an independent procedure. 
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